Scattering theory of multilevel atoms interacting with arbitrary radiation fields 
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We present a generic transfer matrix approach for the description of the interaction of atoms 
possessing multiple ground state and excited state sublevels with light fields. This model allows us 
to treat multi-level atoms as classical scatterers in light fields modified by, in principle, arbitrarily 
complex optical components such as mirrors, resonators, dispersive or dichroic elements, or filters. 
We verify our formalism for two prototypical sub-Doppler cooling mechanisms and show that it 
agrees with the standard literature. 
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I. INTRODUCTION 

The two-level model of atoms interacting with light 
fields [l| has often been used to explore optical cool- 
ing mechanisms @, S 0]- ^s inherent simplicity — the 
atom has one ground state and one excited state — makes 
the resulting models amenable to analysis, but also sup- 
presses mechanisms @ that, in the appropriate parame- 
ter regimes, dominate the interaction. 

A notable example of such an initially overlooked 
mechanism in atomic physics is three-dimensional op- 
tical molasses [(J. By means of the two- level model, 
one can predict the equilibrium temperature, the so- 
called "Doppler" temperature T D , of atoms in molasses 
to be HT, where T is the (half-width at half-maximum) 
linewidth of the transition from the excited to the ground 
level [H . Data from early three-dimensional molasses ex- 
periments contradicted this [7j , showing that the achiev- 
able equilibrium temperature was in fact much lower. 
This discrepancy was resolved independently by two 
groups [f| Q, both explanations relying on the inclu- 
sion of the manifold of magnetic sublevels in each of the 
ground and excited states. In particular, the motion of 
the atoms in the optical field leads to a non-adiabatic fol- 
lowing of the magnetic sublevel populations, which gives 
rise to a strong viscous force and efficient cooling to tem- 
peratures significantly lower than the Doppler tempera- 
ture. 

We recently @ explored a new scattering theory that 
deals with the interaction of light and matter in a uni- 
fied form applicable from microscopic to macroscopic sys- 
tems. In that work we only considered the two-level atom 
model and showed, in particular, how our model can ex- 
plain such mechanisms as standard optical molasses and 
mirror- mediated cooling [loj ] . In this paper we extend 
this model to deal with magnetic sublevels, in much the 
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FIG. 1: Moving scatterer interacting with four field modes 
represented by the Jones vectors A, B, C, and D. The scat- 
terwe has velocity v and is described by means of its polariz- 
ability tensor £. The field mode amplitudes are, in general, 
functions of the wavenumber k. 



same spirit as ref. In due course, this extension will 
enable us to deal with multilevel atoms interacting with 
an arbitrarily complex system composed of immobile mir- 
rors, cavities, MEMS devices, etc., without resorting to 
a quantized model for such a system. 
After we introduce the general extension in the next 
section, we then proceed to explore two prototypical 
systems — the J = \ — ► J' = I transition, leading to the 
"Sisyphus" cooling mechanism, and the J = 1 — > J' = 2 
transition — in Section [Hi] and Section [iVl respectively. 



II. A TRANSFER MATRIX RELATING JONES 
VECTORS 

We investigate the interaction of atoms with light of 
different polarizations. To this end, we denote the two 
polarization basis vectors by fj, and v, whereby the stan- 
dard circular polarization basis is equivalent to setting 
fx = a + and v = o~ . Starting from the transfer ma- 
trix model explored in ref. Q and using the definitions 
in Fig. [TJ we replace each of the field modes by a cor- 
responding Jones vector, similar to the model used in 
ref. Thus, for example, 



A(k) -> A(k) 



A u {k) 



(1) 
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and similarly for B, C and D. The transfer matrix M, 
describing the effect of the scatterer on the four field 
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modes by means of the relation 

(A(k)\ (C(k) 

is now transformed into an order 4 tensor of the form 



(2) 



M 



mn rn 12 
m 21 m 22 



(3) 



where each of m a p (a, (3 = 1, 2) is a 2 x 2 matrix relat- 
ing the respective Jones vector components. A general 
recipe for transforming the formulae for the field mode 
amplitudes, as given in ref. H], can be summarized by 
means of the two replacements 



1 
1 



and £ — > C • 



wherever necessary. In particular, then, 



M 



1-tC -i( 



M = 



i-iC 



F *C 



(4) 



(5) 



We follow ref. [12], Complement Em §3-b, in defining 
the polarizability tensor £ as the steady-state expectation 
value of the polarizability operator x; £ is therefore given 
by the trace 



C = Tr(p st .x)=£<jVWlxl.7}, 



(6) 



where p st is the steady-state density matrix describing 
the system and the summation runs over all the internal 
sublevels of the atom, and where we construct the order 
4 polarizability operator tensor x similarly to ref. [l3|, 
Eq. (14.9-24). In the general /i, v basis: 




i\d v \e) . 



0K|e) 
OKIe) 



(z|d M |e)(e|d Al |j) (i\ | e) (e| d u \ j ) 
{i\d v \e){e\d^\j) {i\d v \e){e\d v \j) 



(7) 



with Co being the characteristic polarizability of the 
atom. In Eq. ([7]), the dipole moment operator d^ (d v ) 
is related to the p {v) polarized light field and the sum 
runs over all the internal sublevels, e, of the atom. The 
matrix elements of d^ (d v ) are given by the appropriate 
Clebsch-Gordan coefficients. 

Importantly, this new transfer matrix still retains all 
its properties, allowing us to model the interaction of the 
multilevel atom with an arbitrary system of immobile op- 
tical elements such as mirrors, cavities, waveplates, etc. 
As in our previous work [9j , this interaction is accounted 
for by the multiplication of the various transfer matrices 
of the elements making up the system; this model is, in 
principle, applicable to systems of arbitrary complexity. 

Finally, we recall that the diagonal elements, (i\p st \i), 
of p st are the populations in each of the sublevels, 




J' = 



FIG. 2: Clebsch-Gordan coefficients for a J 
transition. 



J' 



whereas its off-diagonal elements, (i\p st \j), are the 
respective coherences. The matrix elements of p st are 
obtained from the appropriate optical Bloch equations 
(see, for example, the procedure outlined in ref. [14j]). 
We note here that, through its dependence on p st , M 
depends on the fields that it helps to determine, and 
thus Eq. will in general become a set of nonlinear 
equations. In cases, like the ones considered in the 
following sections, where only one multilevel atom is 
interacting with a linear optical system, this problem 
may be solved using a procedure similar to the one 
outlined below: the fields surrounding the atom are 
obtained from the input fields through linear operations 
and then used with the optical Bloch equations to obtain 
the populations and coherences of the atom's various 
levels. Knowledge of these quantities then determines 
the fields, and hence the forces acting on the atom, 
completely. 

In the following sections we will restrict our discus- 
sion to the case where the input field is not modified by 
other transfer matrices. We will apply this mechanism 
to investigate the behaviour of atoms in two cases where 
the polarization of the light varies in space on scales of 
the order of the wavelength to verify the validity of the 
model given by Eq. to Eq. 0. In the first instance, 
we illuminate our atom with two counterpropagating lin- 
early polarized beams. We choose the planes of polar- 
ization of the two beams to be orthogonal to each other. 
The second configuration we will investigate involves illu- 
minating the atom with two circularly polarized beams, 
choosing opposite handedness for the two beams. These 
two cases mirror those in ref. Q. 



III. ATOMS IN A GRADIENT OF 
POLARIZATION 



In this and the following sections, we will adopt the 
low-intensity hypothesis. This allows us to simplify the 
optical Bloch equations and resulting system consider- 
ably by neglecting the populations and coherences of the 
excited state sublevels. We can thus replace p st by the 
ground state steady-state density matrix, p|*. We de- 
note the diagonal element of p|* , the population in 
sublevel i, by IT, and the off-diagonal element (z, j), the 
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coherence between sublevels i and j, by C,j. 

Here we will discuss what is perhaps the simplest tran- 
sition between two levels with multiple magnetic sub- 
levels: the J = h — ► J' = § transition. In this case, we 
have two ground sublevels so that p S g is a 2 x 2 matrix. 
Fig. [2] tabulates the Clebsch-Gordan coefficients required 
to evaluate £. We thus have: 



2 2 ,+ 2 

c + i_i n i 

• ^2 ' 2 ^2 



(8) 



and 



Co 



whereby 
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(9) 




(10) 



Suppose, now, that we illuminate the atom with two 
counterpropagating beams having orthogonal linear po- 
larization and equal intensity. This can be represented 
by setting 



B(fc) = ^ cxp(ikx - in/A) 



and 



C(fc) 



iB 



■i 



exp(— ikx + z7r/4) 



(11) 



(12) 



where the shift in the x coordinate is introduced to sim- 
plify our expressions. Using the optical Bloch equations, 
we can show that the steady state populations in the 
ground sublevels at zero atomic velocity are given by 




J' = 2 



mj = -1 



FIG. 3: Clebsch-Gordan coefficients for a J — 1 — >J' = 2 
transition. 

where we have assumed that \\kdC/dk\\ ^> ||£||. The 
velocity-dependent force terms in the above expression 
arise through the Doppler shifting of photons both be- 
tween field modes in the same polarization and between 
field modes in different polarizations; these mechanisms 
are accounted for by the diagonal and off-diagonal terms 
in C, respectively. These terms emerge through the 
velocity-dependent terms in the generalised transfer ma- 
trix. 

In the present case, Eq. I|14[) simplifies approximately to 



F = ±hkC,o\B\ 2 sin(2fcx)(n_ ) 



n_i 

2 



assuming that Co is real for simplicity. 
We now let t v be a characteristic residence time of the 
two ground state sublevels; this will introduce a non- 
adiabatic following term, proportional to u, in the pop- 
ulations of each of the sublevels and emerges from the 
optical Bloch equations. Thus, we obtain the expression 



F 



lhk\B\ 2 ( Q sm(Akx) - ^hk 2 \B\ 2 ( vT p sin 2 (2kx), 

(15) 

which agrees precisely with the standard literature (cf. 
Eqs. (4.20) and (4.23) in ref. [|). 



IV. ATOMS IN A GRADIENT OF ELLIPTICITY 



II i = cos (kx) and II i = sin (kx) , 

~2 + 2 



(13) 



noting that the populations do not depend on the field 
amplitudes in the low intensity regime. 

We work to lowest order in Co an d make use of the 
above relations to find the net force acting on the atom: 
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F = hk(\A\ 2 

= 2fifcIm|[C(B + C)] 



C| 2 -|D| 2 ) 
B-C)*} 



4f hk Im{ (CB) • C* + (CC) • B*} 



- 2^hk 2 Im 



— (B 

dk K 



(B + C) 



2fifc|S| 2 Im|[C(B + C)] • (B-C)*} 
2^hk 2 lml\^-(B + C)} ■ (B + C)* 



dk 



(14) 



If we illuminate an atom with two counterpropagating 
beams of light in a <7 + -er~ configuration, rich dynamics 



are obtained not in the simplest (J 



J' = |) 



case, but in the next simplest, where the ground state 
has three magnetic sublevels (J = 1) and the excited 



state five (J' 
and x as 



2). In this case, then, we can express p s q 



P at 



and 



X = Co 



n_i c_io c_i + i 
c+i -i c+i,o n + i 
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using the Clebsch-Gordan coefficients in Fig.[3j Together, 
these give 



C = Co i 



h 

1 



+ 



if 







c 





\ 



n n 



c* 



1 



n 



(18) 



ref. By observing the natural correspondence be- 
tween Co and s± in this latter reference, we can see that 
our expression for the force acting on the atom again 
agrees with the standard literature to first order in | (cf. 
Eq. (5.9) in ref. [5j). The resulting friction force is thus 
due to both the Doppler shift, as evident in the terms 
shown explicitly in Eq. (|2Tj) . as well as through the non- 
adiabatic following of the atomic sublevel populations. 



with C = C + i,-i = C*_ 1+1 = (+l|Pg |— 1) representing 
the nonzero coherence between the mj = +1 and the 
rn / = — 1 sublevels. Note that we again apply the low 
intensity hypothesis, thereby replacing p st with p*£. 

We now illuminate the atom with two counterpropa- 
gating beams of equal intensity, B and C, possessing cr + 
and a" polarization, respectively: 



B(fc) = B[h exp(ikx) 



and 



C(fc) = B ( " ) exp{-ikx) . 



(19) 



(20) 



We again use Eq. (Tl4^) to derive the force acting on the 
atom, which is given by 

F = 2tt|B| 2 Im{|Co(n +1 -n_i) 

+ li( lm{C exp(-2ikx)}} 

- 2^hk 2 \B\ 2 imid^/dk] (j(n +1 + n_i) + n 

+ ±Re{Cexp(-2ifcr)} 
(21) 

where the populations and coherences are again obtained 
from the optical Bloch equations, and can be found in 



V. CONCLUSIONS 

By revisiting the transfer matrix formalism and ex- 
pressing the polarizability of a scatterer as the expec- 
tation value of a quantum operator, we have endowed it 
with a strong quantum character that allows us to handle 
atoms with multiple ground and excited state sublevels. 
In principle, our extended formalism is only limited by 
its reliance on the optical Bloch equations to give expres- 
sions for the ground state populations and coherences; we 
have retained the character of our earlier formalism that 
allowed us to work to arbitrary order in the polarizabil- 
ity. We have applied this theory to two standard sub- 
Dopplcr cooling configurations, the so-called "lin-_L-lin" 
and "<7 + -<7 - " configurations, and thereby reproduced the 
known expressions for the force acting on the atom. 
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